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Polycomb group (PcG) proteins are essential regulators of hematopoietic stem cells. Recent 
extensive mutation analyses of the myeloid malignancies have revealed that inactivating 
somatic mutations in PcG genes such as EZH2 and ASXL1 occur frequently in patients with 
myelodysplastic disorders including myelodysplastic syndromes (MDSs) and MDS/myelopro- 
liferative neoplasm (MPN) overlap disorders (MDS/MPN). In our patient cohort, EZH2 
mutations were also found and often coincided with tet methylcytosine dioxygenase 2 
{,TET2) mutations. Consistent with these findings, deletion of Ezh2 alone was enough to 
induce MDS/MPN-like diseases in mice. Furthermore, concurrent depletion of Ezh2 and 
Tet2 established more advanced myelodysplasia and markedly accelerated the development 
of myelodysplastic disorders including both MDS and MDS/MPN. Comprehensive genome- 
wide analyses in hematopoietic progenitor cells revealed that upon deletion of Ezh2, key 
developmental regulator genes were kept transcriptionally repressed, suggesting compensa- 
tion by Ezhi, whereas a cohort of oncogenic direct and indirect polycomb targets became 
derepressed. Our findings provide the first evidence of the tumor suppressor function of 
EZH2 in myeloid malignancies and highlight the cooperative effect of concurrent gene 
mutations in the pathogenesis of myelodysplastic disorders. 



Recent genomic studies have identified a series 
of recurrent somatic mutations in patients with 
myeloid mahgnancies, including myelodysplas- 
tic syndrome (MDS), myeloproliferative neo- 
plasm (MPN), chronic myelomonocytic leukemia 
(CMML) — classified as an MDS/MPN— and 
acute myeloid leukemia (AML). Notably, these 
mutations frequently occur in epigenetic regu- 
lator genes such as tet methylcytosine dioxygenase 
2 (TET2), isocitmte dehydrogenase 1 {IDHf), IDH2, 
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DNA methyhransferase 3 A {DNMT3A), and poly- 
comb group (PcG) genes (Chung et al., 2012; 
Raza and GaliH, 2012; Shih et al, 2012). PcG 
proteins fomi the polycomb repressive complexes 
(PRCs) 1 and 2. PRC2 contains three core sub- 
units: SUZ12, one of the BED isofonns, and the 
histone methyltransferase EZHI or EZH2, which 
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Figure 1. Distribution of mutations in epigenetic reg- 
ulator genes in patients with myelodysplastic disorders. 

(A) Mutations of TET2, IDHl/2, EZH2, and ASXLl'm 52 sam- 
ples from 119 patients witin MDS, CMML, and AML/MRC 
sinown by colored bars. Each column represents 1 of the 52 
samples with a mutation in one or more of the genes listed. 
EZH2 mutations include uniparental disomy (UPD) of EZH 
mutation and deletion of EZH2 (located at 7q36) associated 
with -7 and 7q- chromosomal abnormalities. The disease 
types of these 52 patients are also shown. Refractory anemia 
with excess blasts (RAEB) includes RAEB-1 and RAEB-2. 
CMML includes CMML-1 and CMML-2. RCUD and RCMD 
denote Refractory cytopenia with unilineage or multilineage 
dysplasia, respectively. Other diseases include 5q- syndrome, 
Refractory anemia with ring sideroblasts (RARS), and MDS, 
unclassifiable (MDS-U). (B) Venn diagrams showing the 
number of patients who have both TET2ar\d PcG mutations. 
The overlap between mutations in 7f72 and EZH2 is statisti- 
cally significant (P < 0.001). 



EZH2 



n=119 



catalyzes di- and tri-methylation of histone H3 at lysine 27 
(H3K27nie3). Canonical PRCl contains four core subunits, 
PCGF (which can include BMIl or MEL18), CBX, PHC, 
and RINGIA or RINGIB, which nionoubiquitylates his- 
tone H2A at lysine 119 (H2AK119ubl). PRCl and PRC2 
cooperatively function as transcriptional repressors via estab- 
lishing these repressive histone modifications. PcG genes have 
been implicated in the maintenance of self-renewing somatic 
and cancer stem cells (Konuma et al., 2010; Sauvageau and 
Sauvageau, 2010) and have typically been characterized as onco- 
genes, as exemplified by somatic activating mutations of EZH2 
in follicular and diffuse large B cell lymphomas (Morin et al., 
2010). However, EZH2, which is located at chromosome 7q36, 



is frequently involved in chromosomal abnormalities such as 
-7 and 7q-, and inactivating mutations of EZH2 have also 
been identified in patients with MDS, MPN, and CMML — 
all clonal myeloid disorders originating from hematopoietic 
stem cells (HSCs; Ernst et al, 2010; Nikoloski et al, 2010). 
Of interest, other components of PRC2, EED, and SUZ12 
appeared to be mutated in a manner similar to EZH2, although 
the frequencies of their mutations are much lower than those 
o{EZH2 mutations. The PRC2-related gene ASXLl, which 
plays an important role in the recruitment and/or stabihty of 
PRC2 (Abdel-Wahab et al, 2012), has also been shown to 
carry inactivating mutations in patients with MDS, MPN, 
and CMML (Abdel-Wahab et al, 2011; Bejar et al., 2011). 
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Figure 2. Generation of mice reconstituted 
with ref2KW, £z/;2^'^, and Tet2^°l'^'>Ezh2^l^ 
hematopoietic cells. (A) A schematic diagram 
of the experimental process. 10^ E14.5 fetal 
liver cells from Cre-ERT^. Cre-ERT:Tet2'^''l'^°, 
Cre-ERT;Ezh^'^^\ and Cre-ERT:Tet2^°l>^°Ezh2"l« 
CD45.2 mice were transplanted into lethally 
irradiated CD45.1 recipient mice. For deletion of 
Ezh2, 100 pi tamoxifen (10 mg/ml) was intra- 
peritoneally injected once a day for 5 con- 
secutive days at 4 wk after transplantation. 
(B) Levels of H3K27me3 in LSK cells and GMPs. 
LSK cells and GMPs from BM of WT, ret2'<W=', 
Ezh2'^l'', and Tet2^°^'^°Ezh2'^l'^ mice were 
analyzed by Western blotting using an anti- 
H3K27me3 antibody at 4 mo after transplanta- 
tion. Levels of H3K27me3 were normalized to 
the amount of H3 and are indicated relative to 
WT control values. The levels of H3K27me3 in 
WT cells were arbitrarily set to 1. The represen- 
tative data from two independent experiments 
are presented. 
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Although these mutations suggest a tumor suppressor function 
of PRC2-related genes in these diseases, both the impact of each 
PcG mutation and its interplay with coinciding mutations re- 
main largely unknown. 

Loss-of-ftinction mouse models of Tet2, DnmtJa, and Bapl, 
as well as gain-of-function mouse models of IDH1, have been 
reported to recapitulate some aspects of malignancies (Moran- 
Crusio et al, 2011; Challen et al., 2012; Dey et al., 2012; Sasaki 
et al., 2012). However, the role of PcG-related mutants has not 
yet been tested in mice. To understand the contribution of in- 
activating PcG mutations to the development of myeloid ma- 
lignancies, we examined the Ezh2-deficient mice over longer 
periods of time and also tested the impact of concurrent deple- 
tion of Ezh2 and Tet2 on hematopoiesis. 

RESULTS 

Somatic mutations of PcG genes in myeloid dysplasia 

In our cohort of 119 patients with myelodysplastic disorders, 
which includes MDS, CMML, and AML with myelodysplasia- 
related changes (AML/MRC), inactivating mutations in 
EZH2 and ASXLl were detected in 8.4 and 16.8% of patients, 
respectively. Moreover, 3.4% of patients had deletion of EZH2 
(located at 7q36) associated with -7 and 7q- chromosomal 
abnormalities (Fig. 1 A and Table SI). Notably, 57.1% of these 
EZH2 mutations coexisted with TET2 mutations. Conversely, 
34.8% of patients with TET2 mutations had coexisting EZH2 
mutations (Fig. 1 B). These findings suggest a link between 
EZH2 and TET2 mutations in the pathogenesis of myelo- 
dysplastic disorders. 

Deletion of Ezh2 results in enhanced repopulating capacity 
of HSCs and promotes myeloid-biased repopulation 

To decipher the pathological role of inactivating EZH2 mu- 
tations and concurrent inactivation of EZH2 and TET2 genes 
in malignant stem cell disorders, we crossed Cre-ERT;Ezh2^'^ 
mice (Mochizuki-Kashio et al., 2011) and Tet2 gene trap mice 
{Tet2™fKD. shide et al, 2012). Tet2™^™ mice, in which the gene 
trap vector is inserted into exon 2 of Tef2 just before the first 
coding exon, express Tet2 niRNA at levels '^20% of those of 
the WT mice and frequendy die by postnatal day 3 (Shide et al., 
2012). Considering the early death of Tet2^'™ mice and a 
necessity to exclude the influence of the loss of Tet2 and Ezh2 
in BM niche cells, we transplanted E14.5 fetal liver cells from 
Cre-ERT control (WT), Cre-ERT;Tet2^'^, Cre-ERT; 
Ezh^'^, and Cre-ERT;Te!2™'^Ezh2!^'^ CD45.2 mice into 
lethaUy irradiated CD45.1 recipient mice and deleted Ezh2 by 
intraperitoneal injection of tamoxifen at 4 wk after transplan- 
tation (Fig. 2 A). We hereafter refer to the recipient mice re- 
constituted with Tef2KD^KD, Ezh2'^'^, and Tet2^''^Ezh2'^'^ 
ceUs as Tet2™'™, Ezh2'^"^, and Tet2^'™Ezh2'^"^ mice, re- 
spectively. EZH2 is a catalytic component of PRC2 that cata- 
lyzes the methylation of H3K27. As expected, the levels of 
H3K27me3 were markedly reduced upon deletion of Ezh2 
in both Lineage^Sca-l"'"c-Kit''" (LSK) cells, which include HSCs 
and multipotent progenitor cells (MPPs), and in granulocyte- 
macrophage progenitors (GMPs; Fig. 2 B).TET2, a methylcytosine 
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Figure 3. Enhanced repopulating capacity of ret2'^°"*'', Ezh2^l^, 
and Tet2^°"^''Ezh2''l^ HSCs. (A) Competitive repopulating assays. 600 LSK 
cells from WT, Tet2^"!>^°, Ezhl'^l^, and Tet^^^'^^Ezhl'^'^ mice at 12 wk after 
deletion of Ezh2\Nere injected into lethally irradiated CD45.1 recipients 
with 4x10^ CD45.1 competitor BM cells. The percentages of donor-derived 
total CD45.2+ hematopoietic cells in the PB (left) and the percentages of 
total CD45.2+ hematopoietic cells and myeloid (Gr-1+ and/or Mac-I*) and 
lymphoid (B220+, CD4+, or CD8+) cells in the PB at 12 wk after transplanta- 
tion (right) are shown as mean ± SEM (n = 4-6). These data are represen- 
tative data from two independent experiments. (B) Representative flow 
cytometric profiles of donor-derived CD45.2+ myeloid (top) and lymphoid 
cells (bottom) in the PB at 12 wk after transplantation. (C) Percentages of 
donor-derived CD45.2+ hematopoietic cells in the BM LSK HSC/MPP and 
myeloid progenitor fractions. The data are shown as mean ± SEM [n = 4-6). 
(D) Secondary transplantation assays. Total BM cells (2 x 10^) from primary 
recipient mice at 12 wk after transplantation were transplanted into lethally 
irradiated recipient mice without competitor cells. The percentages of 
donor-derived total CD45.2+ hematopoietic cells and myeloid and lym- 
phoid cells in the PB at 20 wk after transplantation are shown as mean ± 
SEM (n = 3-5). *, P < 0.05; **, P < 0.01 ; ***, P < 0.001. 



dioxygenase, catalyzes the oxidation of 5-mC (5-methylcyto- 
sine) to 5-hmC (5-hydroxymethyl cytosine), the first step of 
active demethylation ((Ernst et al., 2010; Ko et al., 2010). The 
levels of 5-hmC in total BM ceUs were also reduced in Tet2^'^ 
and Tet2™'^^Ezh2^'^ cells, whereas there was no obvious 
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Figure 4. Hematopoiesis in the recipient mice 
reconstituted with ref2^W, £z/;2^/^, and Tet:^l^° 
Ezh2^^ hematopoietic cells. (A) PB cell counts in 
recipients reconstituted with ret2'^W, Ezh2^^^, and 
Jet2^°<'^°Ezh2'^"^ fetal liver hematopoietic cells at 12 wk 
after deletion of Ezh2. White blood cell (WBC), hemo- 
globin (Hb), and platelet (Pit) counts in the PB from WT 
(n = 19), ret2KW<D = 23), Ezhl"!^ [n = 13), and 
7gf2KWKD£^^2'^;A ^ 23) mice are plotted as dots and 
the mean values are indicated as bars. (B) The hemato- 
poietic compartment of the PB assessed by manual 
differential cell counts, performed at 12 wk after dele- 
tion of Ezh2. WT (n = 5), Tet^°l«> (n = 7), Ezh2'^l'^ (n = 5), 
and Tet^^^'^^Ezh^l^ [n = 6] mice are represented 
as mean ± SEM. Arrowheads indicate statistically sig- 
nificant populations compared with WT (P < 0.05). 
(C) Smear preparation of PB from WT, Tet^°''^°, Ezh2^^, 
and Tet^*°Ezh^l^ mice at 1 2 wk after deletion of 

observed after May-Giemsa staining. Bar, 10 pm. 
Hyposegmented neutrophils consistent with a pseudo 
Pelger-Huet anomaly (bottom right), hypersegmented 
neutrophils (bottom left), and dysplastic monocyte 
(top right) are depicted in insets. Bar, 5 pm. (D) BM 
analysis of WT, Tet2^°l'^°, Ezh2'^'^,mA Tet2^°l'^"Ezh2^l^ 
mice at 12 wk after deletion of Ezh2. Absolute num- 
bers of total BM cells from a unilateral pair of the 
femur and the tibia of WT (n = 7), Tet2^"!'^'> [n = 9), 

[n = 8), and Tet2^°l'^°Ezh2^l'^ mice (n = 7) are 
plotted as dots and mean values are indicated as bars 
(left). The percentage of LSK cells (middle) and myeloid 
progenitors (right) are presented as mean ± SEM (WT, 
n = 7; ref2^WKD^ „ = g; Ezh2''l^ n = 6; M2^°''^°Ezh2'^^^, 
n = 6). (E) Extramedullary hematopoiesis in the spleen 
of Tet2!^°l'^°, Ezh2^\ and Tet^^l'^^Ezh^" mice at 12 wk 
after deletion of Ezh2. Spleen weights of WT (n = 8), 

mice [n = 7) are plotted as dots and mean values are 
indicated as bars (left). Percentage of LSK cells (middle) 
and myeloid progenitors (right) are presented as 
mean ± SEM (WT n = 5; Tet2^l^, n = S; Ezh2'^\ 
n = 7; ref2^W™fz/72^'A, n = A].*,?< 0.05; 
**, P< 0.01;*", P< 0.001. 
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difference in the levels of 5-niC among each genotype (un- 
published data). 

To explore the consequence of loss of Ezh2 and/or Tet2 
in hematopoietic stem/progenitor cells, we first performed 
competitive repopulating assays using LSK cells recovered from 
the recipient mice at 3 mo after deletion o( Ezh2. It has been 
reported that Tet2 loss leads to increased HSC self-renewal 
and myeloid transformation (Moran-Crusio et al., 2011; Shide 
et al., 2012). As expected, Tet2^^^ LSK cells readily estab- 
Ushed progressively increasing chimerism in the PB over time 
(Fig. 3 A). Similarly, the contribution oiEzhl'^'^ and Tet2^'™ 
Ezli2^'^ LSK-derived cells to the myeloid compartment of the 
PB increased over time, although the total chimerism was lower 
than that of Tef2™^'™ cells because of inefficient production 
of lymphocytes by Ez/!2-deficient LSK cells (Fig. 3, A and B). 
Indeed, LSK cells of all genotypes established significantly 

2630 



higher chimerism in the BM LSK fraction compared with WT 
LSK cells (Fig. 3 C). All mutant cells taken from the primary 
recipients continued to propagate in the secondary recipients 
during serial transplantation and exhibited myeloid-biased dif- 
ferentiation (Fig. 3 D). These results reveal that the loss of Ezh2 
augments the repopulating capacity of HSCs and activates the 
production of myeloid cells in much the same way as the loss 
of Tet2 does. 

Deletion of Ezh2 causes myeloid dysplasia in mice 

We next analyzed the hematopoiesis in recipient mice reconsti- 
tuted with Tcf2™^™, Ezh2^'^, and Tet2'^''^Ezh2^"' ceUs. PB 
from Ezh2'^"^ and Te!2^'^Ezh2''"^ iiuce at 3 mo after dele- 
tion of Ezh2 showed reduced white blood cell counts due to 
lymphopenia and increased platelet counts. In addition, mild but 
significant anemia was detected in these mice (Fig. 4, A and B) . 
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Figure 5. Ezhl^l^ mice develop myelodyspfastic disorders. 

(A) Kaplan-Meier survival curve of WT (n = 24), Tet2^°l'^° 

[n = 28), Ezh2''l'^ [n = 23), and Tet^^^^^Ezhl'^l'^ [n = 32) mice. 
Some of the mice (WT, n = 1 1 ; Tet2^°^'^°, n = 1 1 ; Ezh2'^^\ 
n = 1 1 ; Tet^°l'^°Ezh2''^'^, n = ^3] were sacrificed for analyses 
in the middle of observation (indicated by short vertical bars 
on the curves). These survival curves were generated from 
two to four independent experiments, n.s., not significant. 

(B) Hematoxylin and Eosin (H&E) staining of paraffin- 
embedded sections of spleen. Bars, 500 pm. (C) Spleen 
weights plotted as dots with the mean values indicated as 
bars (WT, n = 3; ret2^W^", n = 4; and 5z/J2'^'^ n = 3). (D) The 
percentage of LSK cells and myeloid progenitors in the 
spleen presented as mean ± SEM (WT, n = S; Tet2^°l'^°, n = 4; 
and fzM-^"^, n = 3). (E) Absolute numbers of total BM cells 
from a unilateral pair of the femur and the tibia are plotted 
as dots and mean values are indicated as bars. (F) Percentage 
of LSK cells and myeloid progenitors in the BM presented as 
mean ± SEM. The numbers of mice analyzed are WT, n = 3; 
Tet2^WD^ n = 4; and Ezh2^^\ n = 3. (G) Hemoglobin contents 
in PB are plotted as dots and mean values are indicated as 
bars (left). Percentage of Annexin V+ cells in Terl ighigi'CD?!^''?'' 
BM erythroblasts is also shown as mean + SEM from three 
independent experiments. The numbers of mice analyzed are 
WT, n = 3; ret2'<W°, n = 3; and fz/72'^''^, n = 3 (right). *, P < 
0.05: **, P < 0.01 ; ***, P< 0.001. 
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Tef2™^'™ mice showed monocytosis in PB as reported pre- 
viously (Moran-Crusio et al., 2011), whereas Tet2^'^Ezh2^'^ 
mice did not (Fig. 4 B). Cytological analysis oi Ezh2^'^ and 
Tet2™'^Ezh2^' ^ PB revealed morphological abnormaHties 
in myeloid cells, such as delayed maturation of neutrophils, 
hyposegmented neutrophils consistent with a pseudo Pelger- 
Huet anomaly, hypersegmented neutrophils, and dysplasia of 
monocytes (Fig. 4 C). This trend was more evident in Te<2™^™ 
Ezh2^'^ PB compared with Ezh2^'^ PB but was not obvious 
m Tet2^'^ PB (Fig. 4 C). 

BM analysis at 3 mo after deletion o( Ezh2 revealed that 
although total BM cell numbers were mildly increased only 
in Tet2^^^ mice, all mutant mice had a greater proportion 
of LSK ceUs than the WT mice (Fig. 4 D). Similarly, Tet2^^^, 
Ezh2^'^, and Tet2^'™Ezh2^'^ mice had a greater propor- 
tion of common myeloid progenitors (CMPs) and/or MEPs 
(Fig. 4 D). Furthemiore, extrameduUary heniatopoiesis was evi- 
dent in the spleen of all mutant mice as defined by a significant 



increase in the absolute number of LSK cells and myeloid pro- 
genitors, although splenomegaly was evident only in Tet2^'™ 
mice, likely because of impaired lymphopoiesis in the absence 
of Ezh2 (Fig. 4 E). 

Ezh2^l^ and Tet2^°l^°Ezh2'^l^ mice develop 
myelodysplastic disorders 

During a long observation period, we found that not only 
Tef2™^™ mice but also Ezh2^'^ mice die of hematological 
disorders. This is the first study to describe the lethal hemato- 
logic phenotypes caused by loss of Ezh2. Tef2'™^™ mice devel- 
oped CMML-like disease accompanied by monocytosis and 
extrameduUary hematopoiesis as reported previously (Moran- 
Crusio et al., 2011) and half of the mice died by 10 mo after 
transplantation. Notably, Ezh2^^^ mice also developed MDS/ 
MPN-like disease and half of the mice died by 10 mo after 
transplantation (Fig. 5 A and Table S2). They showed myelo- 
proliferative features characterized by active extrameduUary 
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Figure 6. Tet^°"^°Ezh2^l^ mice develop myelo- 
dysplastic disorders. (A-C) Hematological data from 
moribund Tet2'^''^'^''Ezh2^^'^ mice and those surviving 
at 300 d after transplantation. PB cell counts and the 
number of monocytes in the PB (A), spleen weights (B), 
and absolute numbers of total BM cells from a uni- 
lateral pair of the femur and the tibia (C, left) are 
plotted as dots and mean values are indicated as 
bars. Percentage of LSK cells (C, middle) and myeloid 
progenitor cells (C, right) are presented as mean ± 
SEM. Data were collected from WT [n = 8), MDS/MPN 
{n = 8), and MDS [n = 6) mice for spleen weights and 
WT (f7 = 1 1], MDS/MPN (n = 7], and MDS [n = 6) mice 
for others. (D) Percentage of Annexin V+ cells in 
Ter119''*CD7l'''9'' BM erythroblasts shown as 
mean + SEM from three independent experiments (WT, 
n = 3; MDS/MPN, n = 2; MDS, n = 2). (E) Percentage 
of Annexin V+ cells in LSK cells from the indicated 
mice at 8 mo after transplantation shown as mean + 
SEM from three independent experiments (WT, n = 3; 
Teti^mo^ n = 3; Ezh2^^'', n = 3; and Tet^^l^Ezhl'^^, 
n = 3). *, P < 0.05; **, P < 0.01 ; ***, P < 0.001. 



hematopoiesis in the spleen as evident from splenomegaly 
with a marked increase in LSK cells and MEPs resulting in ef- 
facement of lymphoid follicles (Fig. 5, B— D). They also showed 
a significant increase in LSK cells and MEPs in the BM, whereas 
they did not show obvious monocytosis in the PB (Fig. 5, 
E and F; and Table S2).They were anemic and showed increased 
apoptosis inTerll9'''e''CD71'"sh erythroblasts (Fig. 5 G) and 
LSK cells (Fig. 6 E) in the BM, suggesting ineffective eryth- 
ropoiesis, a feature compatible with myelodysplastic disorders. 
Ezh2^'^ mice also showed dysplasia of myeloid cells as de- 
scribed above (Fig. 4 C). 

To our surprise, concurrent deletion of Tet2 and Ezh2 sig- 
nificantly shortened the latency of disease development and all 
Tet2^'^Ezh2'^''^ mice died by 10 mo after transplantation 
(Fig. 5 A). Histological analysis of moribund Tet2^^^Ezh2^^'^ 
mice revealed lethal pneumonia in all mice analyzed [n = 6; 
unpublished data). Tet2^^™Ezh2^^^ mice developed hetero- 
geneous diseases which were classified into two types of my- 
eloid neoplasms: (1) MDS/MPN including CMML (« = 8) and 
(2) MDS {n = 6; Table S2). MDS/MPN mice showed myelo- 
proliferative features, including CMML-hke monocytosis in 
the PB (Fig. 6 A and Table S2) and/or splenomegaly with ex- 
trameduUary hematopoiesis (Fig. 6 B), and an increase in LSK 
cells in the BM (Fig. 6 C). In contrast, MDS mice did not 
show obvious myeloproliferative features but showed a trend 
of pancytopenia (Fig. 6 A and Table S2). Myeloid dysplasia, 
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including delayed maturation of neutrophils, a pseudo Pelger- 
Huet anomaly, hypersegmented neutrophils, and dysplasia of 
monocytes, was evident in Tet2}^'^Ezh2^'^ mice compared 
with Tet2^'^ and Ezh2^'^ mice as described above (Fig. 4 C). 
Furthemiore, MDS and MDS/MPN mice had a lower propor- 
tion of Terl 1 9^ erythroblasts compared with WT mice (unpub- 
hshed data). The proportion of Annexin cells in CDTl'"^'' 
Terll9'^s'' erythroblasts was significantly higher in both MDS 
and MDS/MPN mice compared with their WT counterparts, 
implicating enhanced apoptosis as a cause of anemia (Fig. 6 D). 
In contrast, enhanced apoptosis in the absence of Ezh2 was 
canceled by loss of Tet2 in LSK cells at 8 mo after transplanta- 
tion (Fig. 6 E), suggesting differential impact of loss of Ezh2 
and Tet2 on apoptosis in disease-initiating cells. 

Alterations in expression of Ezhi target genes 
associated with disease progression 

To elucidate the changes in gene expression responsible for 
the hematological phenotypes induced by loss of Ezh2 and/ or 
Tet2, we performed niicroarray analysis using LSK cells and 
GMPs from each mutant mouse at 3 mo after deletion o(Ezh2 
and MDS/MPN and MDS mice. Of interest, there was signif- 
icant overlap between genes up-regulated or down-regulated 
in Tet2^'^ and Ezh2'^"^ LSK cells. This was true among all 
genotypes and diseases, and the overall gene expression pro- 
files were highly correlated with each other (Fig. 7). To better 
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Figure 7. Gene expression alterations associated 
with disease progression. (A) Venn diagram of genes 
up-regulated or down-regulated in LSK cells from 
7-et2KD;KD^ Ezh2^l^, and Tet2^l'^°Ezh2^^'^ mice at 4 mo 
after transplantation and MDS/MPN (ID19) and MDS 
(ID18) mice (>2.0- or <0.5-fold compared with the 
WT control, respectively). The numbers of genes in 
each group are indicated. The overlaps between every 
two gene sets are statistically significant (P < 0.01). 
(B) A heat map showing the correlation of overall 
gene expression profiles among LSK cells from the 
indicated mice. Black, red, and blue colors represent 
no, positive, and negative correlation, respectively. 
The score of correlation coefficient between LSK cells 
from each genotype or diseased mice was defined 
with Pearson's correlation and is shown on the right. 



understand the consequences of the loss of Ezh2 and/or 
Tet2, we next performed gene set enrichment analysis (GSEA; 
Subramanian et al., 2005). Embryonic stem (ES) cell— like sig- 
natures, which have been widely used in assessment of cancer 
gene signatures, are subdivided into a PRC module, a Myc mod- 
ule, and a core ES module (Kim et al., 2010). Among these, the 
Myc module, which is active in various cancers and predicts 
cancer outcome, was significantly enriched in Ezli2^'^ LSK cells 
and became highly enriched m Tet2^'™Ezh2'^"^ LSK ceUs 
during the development of MDS/MPN and MDS (Fig. 8 A). 
Of note, the Myc module became enriched in LSK cells but 
not in GMPs, suggesting the existence of pathological stem 
cells in the LSK fi-action. Several modules related to enhanced 
cell cycle were also positively enriched in MDS/MPN and 
MDS LSK cells (unpubUshed data). We next analyzed gene sets 
of Ezh2 targets that lose H3K27nie3 upon deletion of Ezh2 in 
ES cells (ES_Ezh2 targets) and Ezhl targets profiled in Ezh2^'^ 
ES cells (ES_Ezhl targets; Shen et al., 2008). As expected, all 
of the PRC2 gene sets showed a trend of positive enrichment 
(derepression) in Ezh2'^"^ and Tet2^'^Ezh2'^"^ LSK cells and 
GMPs (Fig. 8 A). Notably, however, Ezhl targets were nega- 
tively enriched (repressed) and the degree of positive enrich- 
ment of Ezh2 targets was attenuated in Tet2^'™Ezh2'^'^ LSK 
cells from mice which developed niyelodysplastic disorders. 
These data suggest compensatory action by Ezhl in LSK cells 
during the development of niyelodysplastic disorders. Of inter- 
est, both the Ezhl and Ezh2 targets were positively enriched 
even in Tet2^'^ LSK cells, suggesting some involvement of 
Tet2 in PRC2-mediated gene silencing in HSC/MPPs. 

Alterations in H3K27me3 upon deletion of Ezh2 

To assess the genome- wide effects of loss of Ezh2 and/ or Tet2 
on epigenome, we conducted chromatin immunoprecipita- 
tion sequencing (ChlP-seq) for H3K27me3 in hematopoietic 
cells from WT, Tc!2^'^, Ezh2'''^, and Tet2'^'™Ezh2''"^ 
mice at 4 mo after transplantation. Due to the paucity of LSK 



HSC/MPPs, we used GMPs for ChlP-seq analysis and charac- 
terized the H3K27me3 levels over the region from 5.0 kb up- 
stream to 0.5 kb downstream of transcriptional start sites (TSSs) 
of Reference Sequence (RefSeq) genes (Pruitt et al., 2007). As 
expected, drastic reductions in the levels of H3K27me3 were 
observed in GMPs from Ezh2^'^ and Tet2^'^Ezh2'^'^ mice, 
although considerable levels of H3K27iiie3 remained even in 
the absence of Ezh2 (Fig. 8 B). Among PcG targets (greater 
than fourfold enrichment of H3K27nie3 compared with input 
in WT GMPs) that showed greater than twofold reduction of 
H3K27me3 levels in Tet2™'™Ezh2^'^ GMPs, 121 genes were 
derepressed greater than twofold in both MDS/MPN and MDS 
LSK cells compared with WT. These genes included 24 poten- 
tial oncogenes such as Hmga2 (Oguro et al., 2012), Phx3 (Li 
et al., 2013), and Lmol (Tremblay et al., 2010;Table S3). 

Surprisingly, H3K27iiie3 marks around TSSs became more 
enriched (greater than onefold compared with WT) in a small 
portion of genes upon deletion o£Ezh2, particularly in Ezh2^'^ 
GMPs (Fig. 8, B— D), suggesting a major role for Ezhl on these 
gene promoters. Indeed, these genes overlapped considerably 
with the targets of Ezhl in ES cells (Shen et al., 2008; Fig. 8 E), 
and this trend was more obvious with genes that showed high 
levels of H3K27iiie3 (>35-fold enrichment of H3K27me3 
compared with input in Ezli2^'^ GMPs). More importantly, 
the genes that maintained H3K27iiie3 levels upon deletion of 
Ezh2 were overlapped markedly with genes marked with bi- 
valent histone domains in ES cells (Bernstein et al., 2006; Ku 
et al., 2008; Fig. 8 F). This trend was again more obvious with 
the 124 genes that showed high levels of H3K27me3 (>35- 
fold enrichment of H3K27me3 in Ezh2^"^ GMPs), and most 
of these genes appeared to be canonical developmental regula- 
tors, including genes encoding homeobox, paired-box, T-box, 
forkhead, and Gata family transcription factors and zinc fmger 
DNA-binding proteins (Table S4). Because the amounts of DNA 
mimunoprecipitated from Ezh2^"^ and Tct2^'™ Ezh2^'^ 
GMPs are less than those from WT and Tet2^^^ GMPs, we 
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Figure 8. Alterations in gene expression 
and H3K27me3 levels upon deletion of 
Ezh2. (A) Gene expression alterations associ- 
ated with disease progression. Normalized 
enrichment score (NES) from overall gene 
expression profiles of LSK cells and GMPs 
derived from GSEA is shown as the number 
in each cell. Red and blue colors represent 
positive (up-regulated in the given genotype 
relative to WT) and negative (up-regulated in 
WT relative to the given genotype) enrichment, 
respectively. Concentrated colors show that 
the nominal p-value is <0.05 and false discov- 
ery rate (FDR) is <0.05, which suggests mean- 
ingful enrichment for the given gene sets. Pale 
colors show borderline enrichment. The cells 
which do not meet these criteria show white 
color. (B) Summary of H3K27me3 enrichment 
detected by ChlP-seq analysis. GMPs isolated 
from primary recipient mice at 4 mo after 
transplantation were subjected to ChlP-seq 
analysis using an anti-H3K27me3 antibody. 
The fold enrichment values of H3K27me3 
signals were calculated against the input 
signals (ChlP/input) from 5.0 kb upstream to 
0.5 kb downstream of the TSSs of RefSeq 
genes. The red boxes indicate genes with the 
H3K27me3 enrichment >35-fold in Ezh2^l'^ 
and Tet^°l'^Ezh2^l^ GMPs. (C) Scatter plots 
of H3K27me3 enrichment in Tet^^''^^, Ezh2^l^, 
and Tet^°l'^°Ezh2^l'' GMPs in comparison with 
that in WT GMPs. Gray lines indicate 35-fold 
enrichment, which is the value defined as high 
enrichment in A. The numbers of genes with 
higher and lower fold enrichment than that 
of WT (divided by red lines) are indicated. 
(D) Scatter plots of 1 24 genes enclosed by a 
red square in B showing the H3K27me3 en- 
richment >35-fold relative to WT in fz/j2'^''^ 
and M2^°l'^Ezh2>^l^ (DKO) GMPs. Gray lines 
indicate 1.0, which is the value indicating that 
Ezh2'^l^ and Tet^°l'^°Ezh2'^'^ GMPs showed 
the same fold enrichment as WT. (E and F) The 
proportion of Ezhi target genes (E) and genes 
that are marked with univalent (H3K4me3 or 
H3K27me3) or bivalent histone domains 
(H3K4me3 and H3K27me3; F) identified in ES 
cells were categorized by the changes in 



H3K27me3 levels in Ezh2'^^'^ GMPs. The proportions in all RefSeq genes are shown in the first lane. PRC2 target genes were subdivided into genes with 
reduced H3K27me3 levels (<0.5-fold compared with WT), genes that maintained H3K27me3 levels (>l-fold compared with WT), and genes with the 
H3K27me3 enrichment >35-fold compared with input in Ezh2^^'^ GMPs in B. 



may somewhat overestimate the H3K27me3 levels in Ezh2^'^ 
and Tet2™^™ Ezh2^"^ GMPs. However, quantitative ChIP 
analysis confirmed the trend of ChlP-sequence data (Fig. 9). 
These results suggest that Ezhl plays a crucial role in keeping key 
developmental regulator genes transcriptionally repressed. 

DISCUSSION 

In this study, we first demonstrated that the deletion of Ezli2 
confers a growth advantage to HSCs and results in an MDS/ 
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MPN-like disease in mice, similar to the depletion of Tet.2. 
Ezh2^"^ mice developed a lethal MDS/MPN-Uke disease 
with myeloproliferative features, including enhanced repop- 
ulating capacity of HSCs and extrameduUary hematopoiesis 
in the spleen, and myelodysplastic features, such as anemia 
accompanied by enhanced apoptosis in erythroblasts and dys- 
plasia in myeloid cells. These findings are compatible with 
the recurrent inactivation of EZH2 genes in myelodysplastic 
disorders and strongly support the tumor suppressor fianction 
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Figure 9. The levels of H3K27me3 at promoter regions in GMPs. Visualization of ChlP-seq data of the H3K27me3 levels of several developmental 
regulator genes [Hoxd locus, Gata4, and PaxS), tumor suppressor genes [Cdkn2a and Cdkn2b) and oncogenes (Hmga2 and Pbx3] in GMPs from WT, ref2'^'^'*'°, 
Ezh2'^l'^, and Tet2!^''^'^°Ezh2^l^ mice at 4 mo after transplantation using the Integrative Genomics Viewer (IGV). Schematic diagram of these gene loci indi- 
cates their genomic structures. Exons and untranslated regions are demarcated by large and small black boxes, respectively. Data of ChIP analyses at the 
promoters of selected genes are also depicted. Quantitative ChIP analyses of GMPs from WT and fz/i2'^"^ mice at 6 mo after transplantation were per- 
formed. The relative amounts of immunoprecipitated DNA are depicted as a percentage of input DNA. N.D. indicates not detected. The data are shown as 
the mean ± SEM for triplicate analyses. Regions amplified from the precipitated DNA by site-specific quantitative PGR are indicated by arrows. *, P < 0.05. 



of EZH2. Although no mutations have been identified in 
PRCl genes, we recently reported that concurrent loss of the 
PRCl gene Bmil along with Ink4a/Arf in mice causes patho- 
logical hematopoiesis that mimics human primary myelofi- 
brosis, a disease categorized as an MPN (Oguro et al., 2012). 
Together, inactivation of PcG genes could contribute to the 
pathogenesis of malignant stem cell disorders. 

We further demonstrated that inactivating Ezh2 and 
TET2 mutations coexist in 6.7% of patients with myelodys- 
plastic disorders, and concurrent deletion o{Ezh2 and Tet2 in 
mice accelerates the development of heterogeneous myelo- 
dysplastic disorders including MDS and MDS/MPN. Of note, 



Tet2^^^Ezh2^^^ mice showed more advanced myelodys- 
plastic features compared with Tet2^^^ and Ezh2^'^ mice. 
This is the first study to describe the cooperative effects of the 
two independent mutations in the pathogenesis of myelodys- 
plastic disorders. Extensive mutation analyses of the myeloid 
malignancies have revealed that most MDS and MDS/MPN 
patients have several concurrent genetic mutations, but their 
pathological meanings remain to be addressed. Our findings 
highhght the cooperative effects of the independent muta- 
tions in the pathogenesis of myelodysplastic disorders. It is 
assumed that the deregulated gene expression upon concurrent 
depletion of Tet2 and Ezh2 accounts for enhancement in the 
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myelodysplastic features, such as dysplasia of myeloid cells, 
functional defects in neutrophils, and enhanced apoptosis in 
erythroid cells, although these details require further investi- 
gation. On the contrary, EZH2 mutations are very rare in 
de novo AML and chronic myeloid leukemia in blastic trans- 
formation (CML-BC; Ernst et al, 2010; Chung et al, 2012; 
Raza and Galih, 2012; Shih et al., 2012), and MDS patients 
with EZH2 mutations have much lower risk of transforma- 
tion to AML compared with those with other gene mutations 
(Chung et al., 2012; Raza and GaliU, 2012; Shih et al, 2012). 
hideed, no EzhZ^^'^ and Tet2^^'^Ezh2^^'^ mice developed 
AML in this study. We and other groups have previously re- 
ported that deletion of Ezh2 largely attenuates the leukenio- 
genecity of AML in a mouse model (Neff et al., 2012; Tanaka 
et al., 2012). These findings suggest differential roles of EZH2 
in AML versus MDS and MDS/MPN. 

Although deletion of Ezh2 results in a drastic reduction in 
H3K27me3 levels, considerable levels of H3K27nie3 are re- 
tained. The genes that maintained H3K27me3 levels after loss 
of Ezh2 largely overlap with the targets of Ezhl in ES cells, 
and most of these genes are developmental regulator genes 
that are marked with bivalent histone domains (Bernstein et al., 
2006; Ku et al., 2008). Furthermore, tumor suppressor genes 
such as Ink4a/Arf, a critical PcG target for the maintenance of 
the self-renewal capacity of HSCs (Oguro et al., 2006), and 
PcG-targeted developmental regulator genes, such as Pax5 
(Oguro et al., 2010), were kept closely repressed by reduced 
but significant levels ofH3K27me3 modification (Fig. 9). These 
findings indicate that Ezhl is capable of repressing develop- 
mental regulator genes and tumor suppressor genes, acting to 
maintain pathological stem cells in the absence of Ezh2. Of 
interest, Ezhl was recently reported to maintain HSCs in a 
slow-cycling, undifferentiated state (Hidalgo et al., 2012). To- 
gether, our findings suggest that Ezhl plays a critical role in 
pathological stem cells as well as HSCs. 

In contrast, several PcG targets became derepressed in 
Ezh2'^^'' and Tet2^^^Ezh2''^'' LSK cells. These included 
many direct PcG targets that are potential oncogenes, such as 
Hmga2, Pbx3, and Lmol. Furthermore, the Myc module, 
which is active in various cancers, was positively enriched in 
Ezh2^'^ LSK cells and further enriched m MDS/MPN and 
MDS LSK cells. The derepression of PcG-targeted oncogenes 
in conjunction with the up-regulation of the Myc module 
could function as the major drivers in HSCs for the induction 
of myelodysplastic and myeloproliferative disorders in the ab- 
sence of Ezh2. 

Of interest, Ezh2 knock-in (Ezh2-KI) mice overexpressing 
Ezh2 specifically in hematopoietic cells develop MPN (Herrera- 
Merchan et al., 2012). Ezh2-KI mice show an enhanced HSC 
activity similarly to Ezli2^'^ mice in this study. GSEA showed 
significant positive enrichment of the Myc module in both 
Ezh2'^'^ and Ezh2-KI LSK cells (unpublished data). How- 
ever, Ezh2-KI mice do not show myelodysplastic features. 
Apoptosis in LSK cells is suppressed in Ezh2-KI mice but en- 
hanced in Ezh2^'^ mice. Ezh2-KI LSK cells show enhanced 
colony-forming capacity, whereas Ezh2'^'^ LSK cells show 

2636 



severely compromised proliferation in culture probably due 
to failure in the transcriptional repression of tumor suppressor 
genes such as, pi 9'^'^ (unpublished data), as we have previously 
reported in transformed myeloid progenitors (Tanaka et al., 
2012). These findings clearly discriminate the effects of over- 
expression and loss/hypomorphisni of Ezh2 on heniatopoie- 
sis, even though both results in clonal propagation of HSCs. 

Finally, the functional cross talk between DNA and histone 
modifications is an intriguing issue to understand the epigeiietics 
of myeloid malignancies. Tet2™'^Ezh2^'^ mice developed 
not only MDS/MPN but also MDS, whereas Ezh2'^"^ mice 
developed only MDS/MPN. Moreover, myelodysplastic fea- 
tures were more pronounced in Tet2™'^Ezh2^'^ mice com- 
pared with Ezh2'^"^ mice. These findings suggest that loss of 
Ezh2 and Tet2 cooperate to generate a MDS phenotype in 
mice. However, enhanced apoptosis in the absence of Ezh2 
was significantly attenuated in LSK cells, but not in erythro- 
blasts, upon down-regulation of the Tet2 expression, suggest- 
ing an opposing impact of loss of Ezh2 and Tet2 on apoptosis 
in MDS-initiating hematopoietic stem and progenitor cells. 

The understanding of the pathological significance of each 
somatic gene mutation in myeloid malignancies is just starting 
to be revealed. Nevertheless, our findings clearly indicate that 
mouse models, including our own, are capable of faithfully 
mimicking myelodysplastic disorders and can serve as valuable 
tools for the analysis of epigenomics of myelodysplastic disor- 
ders and preclinical therapeutic studies. 

MATERIALS AND METHODS 

Sequencing analysis of primary samples with MDS and related 
myeloid neoplasms. Genomic DNAs from 119 BM samples with MDS 
and related neoplasms were analyzed for mutations of TET2, IDHl, IDH2, 
EZH2, and ASXLl genes by high- throughput sequencing of SureSelect- 
enriched target exons (Agilent Technologies) according to manufacturer's pro- 
tocol. Genome-wide copy number lesions were also analyzed using GeneChip 
SNP-genotyping microarrays (Aflymetrix) as previously described (Chen 
et al., 2008). Allelic imbalances were detected from the allele-specific copy 
numbers calculated using CNAG/AsCNAR software (http://www. genome 
.umin.jp; Nannya et al., 2005; Yamamoto et al., 2007). AH mutations were 
confirmed by Sanger-sequencing. AH patient samples were appropriately 
anonymized before the analysis according to the protocol approved by the 
Ethics Boards of the University of Tokyo. 

Mice and generation of fetal liver chimeras. Tet2^^^ and Ezh^^^ 
mice, which had been backcrossed at least eight times onto a C57BL/6 
(CD45.2) background, were used (Mochizuki-Kashio et al., 2011; Shide 
et al., 2012). For conditional deletion of Ezh 2, the mice were crossed with 
Rosa: :Cre-ERT mice (Taconic). C57BL/6 (CD45.2) mice were purchased 
from Japan SLC. C57BL/6 mice congenic for the Ly5 locus (CD45.1) were 
purchased from Sankyo-Lab Service. All experiments using the mice were 
performed in accordance with our institutional guidelines for the use of 
laboratory animals and approved by the Review Board for animal experi- 
ments of Chiba University (approval ID: 25—104). To generate fetal liver 
chimeras, we intravenously transplanted E14.5 fetal liver from CD45.2 mu- 
tant mice (test cells) into 8-wk-old CD45.1 recipients irradiated at a dose of 
9.5 Gy. After 1 mo, the recipient mice were injected with 100 pi tamoxifen 
dissolved in corn oil at a concentration of 10 mg/ml intraperitoneally once 
a day for 5 consecutive days to induce Cre activity. 

Serial transplantation assay. 600 CD45.2 LSK cells were sorted from the 
BM of the primary recipient mice reconstituted with WT, Tct2^^'^, Ezh2^^^, 
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and Tet2^^^Ezh2^^^ cells at 4 mo after transplantation and injected into 
lethally irradiated CD45.1 niice (secondary recipient mice) along with 4 x 
10^ BM cells (CD45.1) for competitor cells. For tertiary transplantation, 2 x 
10^ whole BM cells taken from the secondary recipient mice 4 mo after 
transplantation were injected into lethally irradiated CD45.1 mice. 

Western blot analysis of H3K27me3. Samples were separated by SDS-PAGE, 
transferred to a PVDF membrane and detected by Western blotting using the 
following antibodies: anti-H3 (Abeam) and anti-H3K27me3 (MUlipore). 

Flow cytometry and antibodies. niAbs recognizing the following anti- 
gens were used in flow cytometry and cell sorting: CD45.2 (104), CD45.1 
(A20), Gr-1 (RB6-8C5), CDllb/Mac-1 (Ml/70), Ter-119, CD127/IL-7Ra 
(A7R34), B220 (RA3-6B2), CD4 (L3T4), CD8a (53-6.7), CD117/c-Kit 
(2B8), Sca-1 (D7), and CD16/32/Fc7RlI-IIl (93). The niAbs were purchased 
from BD, eBioscience, or BioLegend. For Annexin V staining, cells were 
suspended with 1 X Annexin binding buffer (BD) and stained with FITC- 
Annexin V (BD) according to the manufacturer's protocol. Dead cells were 
eUminated by staining with 1 |J.g/ml propidium iodide (Sigma-Aldrich). All 
flow cytometric analyses and cell sorting were performed on a FACSAria II 
or FACSCanto II (BD). 

Microarray analysis. A one-color microarray-based gene expression analy- 
sis system (Agilent Technologies) containing 39,429 clones was used (Sure- 
Print G3 Mouse GE 8 x 60K array; Agilent Technologies), according to the 
manufacturer's instructions. Total RNA was extracted from 5 X 10''^ sorted 
LSK cells or GMPs. A total of 10 ng of total RNA was mixed with spike-in 
controls using a One Color Spike Mix kit (Agilent Technologies), amplified 
and labeled with Cyanine 3 using a Quick Amp Labeling kit (Agilent Tech- 
nologies) according to the manufacturer's instructions, which generated 
single-color labeled cRNA. A total of 1,650 ng of the labeled cRNA was 
used for each hybridization. The process of hybridization and washing was 
perfomied using a Hi-RPM Gene Expression Hybridization kit (Large; Agilent 
Technologies) and a Gene Expression Wash Pack (Agilent Technologies), 
respectively. A DNA microarray scanner (Agilent Technologies) was used 
for array scanning. 

qRT-PCR. Total RNA was isolated usmg TRIZOL LS solution (Invitro- 
gen) and reverse transcribed by the ThemioScript RT-PCR system (Invitro- 
gen) with an oHgo-dT primer. Real-time quantitative RT-PCR (qRT-PCR) 
was performed with an ABI Prism 7300 Thermal Cycler (Apphed Biosystems) 
using FastStart Universal Probe Master (Roche) and the indicated combina- 
tions of Universal Probe Library (Roche) and primers Hsted below. 

Primer sequences and probe numbers for real-time RT-PCR (forward/ 
reverse): Hmga2, 5'-AAGGCAGCAAAAACAAGAGC-375'-CCGTTTT- 
TCTCCAATGGTCT-3'; Phx3, 5'-GCCTTGGAGCAAACTCACTG-37 
5'-AGATGGAGTTGTTGCGTCCT-3'; Hprtl, 5'-TCCTCCTCAG- 
ACCGCTTTT-375'-CCTGGTTCATCATCGCTAATC-3'. Primer 
sequences for manual ChIP analysis at promoter regions: (forward/reverse 
primer) HoxdlO, 5'-TCTCATTGGCTTGGTTGTCA-375'-CAGGAGA- 
GCTGTTGGGAAAG-3'; HoxdP, 5'-TTAAACAACACGCCAAGCTG-37 
5'-AAACCCATCCCTACACACACA-3'; Gat(i4. 5'-ACAGCTCTCT- 
GGTGGCTCTC-375'-AGGTGGTATTCCAGCCCTCT-3'; Pax5. 
5'-CCACTGGACAATGGCAGTTT-375'-TCCCCCTATCCTCG- 
AACTCT-3'; Cdkn2a, 5'-GATGGAGCCCGGACTACAGAAG-37 
5'-CTGTTTCAACGCCCAGCTCTC-3'; Hm£a2, 5'-AAAACTTGGG- 
CTCCGGGTGCAGA-375'-GGGCGCCCAGCTCAGCTCTAG-3'. 

Chip assay and ChIP assay coupled with massive parallel sequenc- 
ing (ChlP-Seq). GMPs from the BM of recipient mice were cross-hnked 
with 0.5% formaldehyde for 5 min at room temperature, washed three times 
with PBS, and lysed with RIPA buff'er (10 niM Tris, pH 8.0, 140 mM 
NaCl, 1 niM EDTA, 1% Triton X-100, 0.1% SDS, 0.1% sodium deoxycho- 
late [DOC], and protease inhibitor cocktail) and sonicated for 30 min using 
a Microson XL2000 Ultrasonic cell disruptor (Misonix). The mean chromatin 



fragment size after sonication was ^^200— 300 bp. After centrifugation, the 
soluble chromatin fraction was recovered, precleared for 1 h at 4°C with a 
mixture of Protein A— and G— conjugated Dynabeads (Invitrogen) blocked 
with BSA and salmon (Invitrogen), and then incubated with an aiiti-H3K27me3 
antibody (07—449; MiUipore) for 2 h at 4°C. Chromatin was UTununopre- 
cipitated overnight at 4°C with antibody-conjugated Dynabeads. The ini- 
munoprecipitates were extensively washed with the following combination 
of wash buffbrs: high salt RIPA buffer (10 niM Tris, pH 8.0, 500 niM NaCl, 
1 niM EDTA, 1% Triton X-100, 0.1% SDS, 0.1% DOC, and protease in- 
hibitor cocktail), LiCl wash buffer (10 niM Tris-HCl, pH 8.0, 250 niM 
LiCl, 1 niM EDTA, 0.5% NP-40, and 0.5% DOC), and TE buffbr (10 niM 
Tris-HCl, pH 8.0, and 1 niM EDTA. Bound chromatm and mput DNA 
were placed in elution buffer (10 niM Tris-HCl, pH 8.0, 5 mM EDTA, 300 niM 
NaCl, and 0.5% SDS) and reverse cross-linked. Ininiunoprecipitated DNA 
and input DNA were treated with RNase A (Sigma- Aldrich) and proteinase K 
(Roche), and purified with a QIAquick PGR purification kit (QIAGEN). 
For Chip assay, quantitative PGR was performed with an ABI prism 7300 
Thermal Cycler using SYBR Premix Ex Taq II (Takara Bio). The primer 
sequences are listed above. For ChlP-Seq, iiimiunoprecipitated DNA and 
input samples were then prepared using ChlP-Seq Sample Prep kit (lUu- 
niina). Adaptor-ligated DNA fragments were size fractionated by 12% acryl- 
amide gel, and the 170- to 250-bp ffaction was recovered. DNA thus 
obtained was amplified by 18 cycles of PGR. 1 ng DNA was used for the 
sequencing reaction of the GAIIx (Illumina), according to the manufactur- 
er's instructions. A total of 170,000—250,000 clusters were generated per 
tile, and 36 cycles of the sequencing reactions were performed. Short-read 
sequences were aligned to the mouse genome sequences (mm9 from Uni- 
versity of California, Santa Cruz Genome Browser; http://genome.ucsc 
.edu/) using the Eland program. Sequences allowing no more than two 
mismatches per sequence were used for the analysis. The ChlP-Seq signal 
was quantified as total number of reads per million. To enumerate genes 
with H3K27me3 modification, genes with cumulative tag numbers in the 
region 5.0 kb upstream of theTSS to 0.5 kb downstream of theTSS greater 
than twofold was selected. 

Statistical analysis. Heat map was drawn by Matrix2png (Pavlidis and 
Noble, 2003). Statistical tests were performed using Prism (version 5; 
GraphPad Software). For analysis of survival curves, log rank (Mantel-Cox) 
test was performed. 

Deposition of the data. Microarray and ChlP-sequence data were depos- 
ited in Gene Expression Omnibus (accession no. GSE42666) and DNA Data 
Bank of Japan (DDBJ, accession no. DRA000485-488), respectively. 

Online supplemental material. Table SI lists the gene mutations and ab- 
normahties in chromosome 7 in patient samples. Table S2 provides the de- 
tailed hematological characteristics of diseased mice. Table S3 hsts the 
potential oncogenes of direct PcG targets derepressed in the absence of Ezh2. 
Table S4 lists the genes with a fold enrichment of H3K27me3 >35 compared 
with input in £'2'/i2-deficient GMPs. Online supplemental material is avail- 
able at http://www.jem.org/cgi/content/full/jem.20131144/DCl. 
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